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We study to what extent a measurement of the m± spectra for hadrons and their resonances can 
resolve ambiguities in the statistical model description of particle production. We describe in a 
quantitative analysis how physical assumptions about the freeze-out geometry and dynamics influ- 
ence the particle spectra. Considering ratios of m±_ distribution of resonance-particle ratios (such 
as K*/K, T,* / A, rj' /rj) we observe significant sensitivity to fireball freeze-out geometry and flow 
dynamics. 

PACS number(s): 12.38.Mh, 25.75.-q, 24.10.Pa 



cn 
o 
o 

(N 
< 

> 
O 

O 



X 



I. INTRODUCTION 

The Fermi statistical model of particle production 1, 
has been used extensively in the field of relativistic 
heavy ion collisions ||5|j . Particle abundancies and spectra 
both at Super Proton Synchrotron (SPS) 0,111111 El 
and Relativistic Heavy Ion Collider (RHIC) [l^ 

in El El 1^ energies have been analyzed in this way. 
The quality of fits to experimental results was such that 
it became possible to discuss hadronization conditions 
quantitatively, but the conclusions of the groups differ. 
For example the values of temperature ran ge f rom as low 
as 110 MeV H El El E3 to 140 MeVjllJllto as high 
as 160, 170, and 180 MeV H i, El Elllllll at both 
SPS and RHIC energies. 

These differences are on a closer inspection not very 
surprising, since the (tacit) assumptions made about 
hadronization mechanisms differ. However, this means 
that before we can say that the freeze-out temperature 
has been determined, we must understand precisely the 
origins of these differences, and proceed to ascertain 
which model is applicable. We shall suggest experimental 
observables which will be particularly sensitive to the dif- 
ferences between the hadronization scenarios, in the hope 
that further experimental study will allow to understand 
the statistical hadronization mechanism. 

We begin with an overview of the differences between 
hadronization scenarios and their relation to the physi- 
cal assumptions used. Every model discussed here has 
been extensively studied before, and has gained accep- 
tance of some part of the heavy ion community. The 
theoretical principles that we invoke are well understood, 
and the methods we use can be found scattered in liter- 
ature. We shall concentrate here on an analysis of reso- 
nances produced in a heavy ion collision. Direct detection 
of hadronic short-lived resonances has become possible 
through invariant mass reconstruction 0, |2l |^ |2^ . 
Resonances have already been proposed as candidates 
for differentiating between freeze-out models Res- 
onances are a sensitive probe of the freeze-out tempera- 
ture, since the ratio of yields of particles with the same 
quark composition is insensitive to both fugacities and 
phase space occupancies, and mass differences are greater 
than the hadronization temperature considered. 



Here we shall develop this reasoning one step further: 
If particles were to be emitted from a static thermal 
source, and feed down corrections were performed, the 
ratio of resonance to daughter particle would be indepen- 
dent of mj_. That this ratio is in general somewhat m± 
dependent is in this situation due to dynamical effects, 
such as hadronizing matter flow and freeze-out geometry. 
For this reason, it can be expected that the m± depen- 
dence of this ratio can help isolate these effects and thus 
remove the ambiguities in the present freeze-out models. 

In Sec. II we review hadronization models and dis- 
cuss their ambiguities. We then show in Sec. Ill how 
resonance m± ratios can be used to distinguish between 
particle hadronization models. We close with a short dis- 
cussion of open issues. 



II. STATISTICAL HADRONIZATION 



A. General remarks 

Nearly all hadronic spectra comprise a significant 
0(50%) component from resonance decays. Fits to data, 
which are not allowing for the decay contributions have 
a very limited usefulness. Particle spectra and thus 
yields are in general, controlled by the properties of 
statistical hadronization models. However, some recent 
work fits the particle slopes only 0, H El El, treat- 
ing the normalization of each particle as a free parame- 
ter. This approach can be argued for assuming a long- 
lived posthadronization "interacting hadron gas phase" 
in which individual hadron abundances subject to in- 
elastic interactions evolve away from chemical equilib- 
rium. This particular reaction picture clashes with e.g. 
the fact that short-lived resonance ratios can be de- 
scribed within the statistical hadronization model using 
the chemical (statistical hadronization) freeze-out tem- 
perature obtained in stable particle studies '18]. This 
implies that in principle the relative normalization of the 
particle spectra should be derived from a hadronization 
scenario involving flavor chemical potentials. In fact a 
study of RHIC spectra finds that the normalization can 
be accounted for |l3 , and that the chemical equilibration 
temperature also describes particle spectra well. This 
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is suggesting that any posthadronization reinteraction 
phase is short and has minor influence on the particle 
yields. 

The problem is that the different ways to derive 
hadronization particle distributions have a profound ef- 
fect on the resulting fitted temperature. Temperature 
affects the absolute number of particles through several 
mechanisms and anticorrelates with the phase space oc- 
cupancy parameters 7^,1 = w, d, s 0, 0| . It has been 
found that the introduction of these parameters, mo- 
tivated by the need to conserve entropy at hadroniza- 
tion decrease the j per degree of freedom consid- 
erably and lowers the freeze-out temperature by 30 MeV 
[T^j l . Other workers assume the light flavors are in chem- 
ical equilibrium 

Considering ratios of resonances to ground state par- 
ticles eliminates the fitted temperature's sensitivity to 
chemical equilibration, since the numerator and denom- 
inator have the same quark composition. In every 
hadronization model considered here the chemical pa- 
rameters cancel out and only temperature and fireball 
freeze-out geometry and dynamics influence the observed 
ratios. 

When fitting the particle spectra, the system's spa- 
tial shape and the way the freeze-out progresses in time 
have a considerable effect on the form of particle distri- 
butions, and hence on the fitted temperature and matter 
flow. The impact of freeze-out geometry and dynam- 
ics on particle spectra were examined well before RHIC 
data became available ^24, 2^ and it was realized that 
an understanding of freeze-out is essential for the statis- 
tical analysis of the fireball [2^ . Even though this matter 
has been clearly recognized, a systematic analysis of how 
freeze-out geometry affects particle distributions is for 
the first time attempted here. In fact, each of the models 
used in the study of particle spectra 0,[l3>ll3>lll>ll3 em- 
ploys a different choice of freeze-out geometry, based on 
different, often tacitly assumed, hadronization scenarios. 
Thus an understanding for the influence of hadronization 
mechanism is impossible to deduce from this diversity. 

However, every study of firball hadronization we are 
aware of uses the Cooper-Frye formula |23 : 

- / x)e{p^d^j:^), (i) 

where is the particle's four-momentum, is the sys- 
tems velocity profile, T is the temperature, A is a chem- 
ical potential, f{E, T, A) is the statistical distribution of 
the emitted particles in terms of energy and conserved 
quantum numbers and describes the hadronization 
geometry. It is the covariant generalization of a vol- 
ume element of the fireball, i.e., a "3D" surface in space- 
time from which particles are emitted. 6{p^^(PYi^) is the 
step function, which eliminates the possible inward emis- 
sion mm. 

The Cooper-Frye formula, Eq. (Q, is believed to be 
the most general way to implement statistical hadroniza- 
tion emission of particles. For it to represent a physical 



description of the system, the following two conditions 
have to be met: 

(i) Statistical hadronization must apply. The parti- 
cles emitted from a volume element (in it's co- moving 
frame) , will be distributed according to the Bosc-Einstein 
or Fermi-Dirac distributions f{E,T, A) for some temper- 
ature T and fugacity A. 

(ii) A "small" volume element hadronizes rapidly in 
it's rest frame, that is, no long lived mixed Quark gluon 
plasma (QGP) -hadronic confined phase exists. 

If this second condition is satisfied, it becomes pos- 
sible to define a hadronization hypersurface = 
{tf{x,y,z),x,y,z)) which specifies at which time tf 
hadrons are emitted from the point {x,y,z). In this 
fast hadronization case differing can be considered 
for physically differing models. The different choices of 
correspond to physically different scenarios, and it 
becomes possible, in principle, to distinguish them ex- 
perimentally. However, if a long-lived mixed phase does 
exist, it might well be that the Cooper-Frye formula can 
be used as an approximation technique to transform a 
hydrodynamically evolving system into hadrons and au- 
thors who worked with a long mixed phase have chosen 
this approach, see e.g. Ref. [3(]| . 

B. Freeze-out geometry 

The high baryon stopping power observed at SPS en- 
ergies j31||, (321, |33|| has prompted some authors to use a 
spherical expansion and freeze-out as an ansatz |lClj| . 

However, at RHIC collision energies the measured 
dN/ drj 13 [s^ indicates that around mid-rapidity the 
system conditions can be approximated by the Bjorken 
picture [3^ . 

To describe particle spectra measured around midra- 
pidity, therefore, boost invariance becomes the dominant 
symmetry on which freeze-out geometry should be based. 
To construct such a hadronization scenario, we consider 
that the most general cylindrically symmetric flow proflle 

/ cosh(yL) cosh(?/^) \ / m_L cosh(y) \ 

^ I sinh(yj_) cos(6l) = ^"^'^'^^ 

sinh(y_L)sin(6') I ' I P-Lsin(0) I 

V sinh(yL) cosh(?/j^) / \m±_ sinh(y) / 

(2) 

(the last, longitudinal coordinate is defined along the 
beam direction) leads to the following rest energy [38| : 

p^u^ = mj^ cosh(?/_L) cosh(?/-?/i)-sinh(?/_L) cos(6'-(/))pj_. 

(3) 

The requirement for the Bjorken picture is that the emis- 
sion volume element has the same j/l dependence: 

p^d3j^^_^cosh(y-t/L) + S. (4) 

This constrains the freeze-out hypersurface to be of the 
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form 



= {tf cosh(?/L), X, y, tf sinh(?/L)). 



(5) 



Here tf is a parameter invariant under boosts in the z di- 
rection, whose physical significance depends on the model 
considered. 

For central collisions, a further simplifying constraint 
is provided by the cylindrical symmetry, which forces t f , 
as well as ijL and to be independent of the angles 
and (j). The freeze-out hypersurface can be parametrized, 
in this case, as 



{t f (r) cosh(yL ) , r sin(6') , r cos(6') , t/ (r) sinh(?/L)) 



d^S'' = tfrdrdOdvL 



dr 



cosh(2/L)^ cos(6'), ^ sin(6'), sinh(yL; 



dtf_ 
dr 



(6) 



(7) 



And the emission clement takes the form 



m±_ cosh(y - hl) ~ P± cos 



tfrdrdOdyL, (8) 



with the same dependence on the angle as Eq. . Equa- 
tion Q can then be integrated over all the possible val- 
ues of j/L and 9 — (j) to give a particle spectrum depending 
purely on the transverse mass, temperature, and y±. The 
fits in ^14,^15, 

are based on such an ansatz. 

What distinguishes the models currently considered is 
the time component of the freeze-out surface. The most 
general freeze-out hypersurface compatible with cylindri- 
cal symmetry is provided by Eq. ( (SJ . Generally, i/ (a 
generic function of r) represents the time, in a frame co- 
moving with the longitudinal flow, at which the surface 
at distance r freezes out. 

The fits in Refs. Hill El are based on a particular 
case of such a freeze-out surface, in which tf is completely 
independent of r {dtf/dr — 0). Such a picture's physical 
reasonableness can be questioned, e.g., why should spa- 
tially distant volume elements, presumably with different 
densities and moving at different transverse velocities, all 
freeze out simultaneously in a longitudinally comoving 
frame?. However, such a simple model can perhaps serve 
as an approximation. 

More generally, the "burning log" model 113 
(sometimes referred to as "blast wave" ; This term, how- 
ever, is also used to refer to the dtf/dr = model de- 
scribed in the preceding paragraph) assumes that the 
emission occurs through a three-dimensional hadroniza- 
tion surface which is moving at a constant "velocity" 
{1/vf — dtf/dr throughout the fireball. Both boost- 
invariant and spherically symmetric versions of burning 
log model were considered. Even if the hadronization 
velocity encompasses an extra parameter Vf, it is worth 
considering since it is based on a physically motivated 



hadronization picture. Moreover, the burning log picture 
is a suitable framework in the study of sudden hadroniza- 
tion. Sudden hadronization occurs when the fireball en- 
counters a mechanical instability [39l |. which combined 
with the fireball's high transverse flow ensures that the 
emission surface spreads to the interior of the fireball 
with ~ c. All of the indications suggested for such 
a picture seem to be borne out by both SPS and RHIC 
data [Ullilill. 

An approach based on the hypothesis of initial state 
"synchronization" by the primary instant of collision and 
the following independent but equivalent evolution of all 
volume elements assumes that each element of the sys- 
tem undergoes freeze-out at the same proper time r. In 
this framework each fireball element expands and cools 
down independently, hadronizing when it's temperature 
and density reach the critical value. This model was suc- 
cessfully used to describe RHIC mj^-spectra ^3- In this 
approach in Eq. © is equal to rcosh(y_L) and the 
hadronization hypersurface in Eq. Q becomes propor- 
tional to the fiow vector: 



d^E^' = TrdrdOdyLu'' 
r = Tsinh(?/j^). 



dVut" 



(9) 
(10) 
(11) 



In this hadronization model the heavy ion fireball be- 
haves similarly to the expanding Hubble universe. In the 
'Hubble' scenario, the Cooper-Frye formula reduces to 
the Touscheck Covariant Boltzmann distribution |^ IbM 

El El. 



(27r)3 ' 



-E/T 



(27r)3 



-p,«VT(i2) 
(13) 



(Where V is the co moving fireball's volume element in 
the local rest frame.) 

To summarize and illustrate the diversity of distinct 
hadronization geometries we present in TableUand Fig.Q] 
the freeze-out scenarios examined here. As we shall see 
the choice of freeze-out geometry produces in a fit of ex- 
perimental data a non trivial effect capable of altering 
significantly the understanding of statistical hadroniza- 
tion parameters. 



C. Flow profile 

Hydrodynamical expansion of the fireball implies in 
general that each volume element will have a different 
density and transverse expansion rate. For this reason, 
the integral over d^E can span a range of flows, weighted 
by density. In flrst approximation one can flt data us- 
ing just an "average" flow velocity throughout the entire 



Longitudinal distance from center { z=Tsinh(yL) ) 



Transverse radius 



FIG. 1: (Color online) While boost-invariance fixes the longitudinal freeze-out structure (left), several scenarios exist for the 
transverse dependence of freeze-out (right). For spherical freeze-out, only plot on the right applies 



TABLE I: (Color online) Freeze-out hypersurfaces at contours of constant radii. 
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»/? = cosh(y^)/T, o = sinh{yj_)/T 



fireball 0,113: 

" J''dr{E-p^^)f{T,yJr),X) 

cx {E-pAf{T,{y,_),X). (14) 
ar 

However, if one wants to properly identify tlie flow 
profile should be taken into account ^45] . Hydrodynamic 
simulations js^l accompanied by assumption that freeze- 
out happens when a volume element reaches a critical 
energy density indicate that the transverse rapidity will 
depend linearly with the radius i.e. v± ^ tanh(r). This 
condition, however, is appropriate for a static freeze-out 



and will not in general hold if the freeze-out is sudden. 
Other flow profiles have been tried in the literature, aris- 
ing from dynamical hypothesis. For example, the as- 
sumption that the freeze-out occurs at the same time tf 
results in a quadratic {v oc r"^) flow profile which has 
also been used recently in fits to data ^3 • the Hubble 
fireball the freeze-out conditions will also result in 
a distinctive flow profile. Specifically with S'' cx; u^, we 
have 7W oc r. 

Density profiles also depend on the assumed initial con- 
dition and the equation of state of the expanding QGP. 
It has been shown p6] that different density choices have 
a considerable effect on both the temperature and fiow 
fits at SPS energies. 
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FIG. 2: (Color online) n,K,K and H m±_ distributions obtained with different freezeout models and flow profiles. For this and 
subsequent figures, a uniform density profile was assumed 



Figure. |21 shows how the choice in hadronization dy- 
namics and flow profiles at same given freeze-out tem- 
perature and transverse flow can result in a range of 
inverse spectral slopes. Here the density profiles were 
assumed to be uniform. It is clear that the same freeze- 
out parameters give rise to a variety of substantially dif- 
ferent particle spectra. Conversely, fits to experimental 
data will only produce reliable information on the freeze- 
out conditions if and when we have a prior knowledge 
of the hadronization geometry and dynamics. There- 
fore, conclusions about statistical model fits, as well as 
arguments whether freeze-out occurs simultaneously for 
different particles or not, cannot be answered while the 
models used to fit the data are plagued by such uncertain- 
ties. We will now turn to study how the measurement of 
spectra of short-lived resonances might provide us with 
a way of making progress. 



III. MOMENTUM DEPENDENCE OF THE 
RESONANCE-PARTICLE RATIOS AS A 
FREEZE-OUT PROBE 

We have shown that the measurement of resonances 
can probe both the hadronization temperature, and the 
lifetime of the interacting hadron gas phase [13. |23 | . Ra- 
tios of a generic resonance (henceforward called Y*) to 
the light particle (which we will refer to as Y) with an 
identical number of valence quarks are particularly sen- 
sitive probes of freezeout temperature because chemical 
dependence cancels out within the ratio. If we examine 
this ratio within a given > m^* range, we expect to 
disentangle flow and freeze-out conditions, since the ratio 
Y* /Y should not depend on mj_ for a purely static and 
thermal source. 

We therefore take the most general Boost-invariant 



freeze-out hypersurface in the Boltzmann limit (see Ta- 
ble ^ Boost invariant implies this is a good approxima- 
tion at midrapidity) 



dN 
drn^ 



oc S{m±,pj_) ^ / rdrS{m±,p±,r), 



(15) 



where 



dtf 

- --^p^Ko{pm^)hiap^), (16) 
or 



with 



cosh[j/^(r)] sinh[j/j_(r) 
P — — , Oi 



T 



T 



(17) 



and use it to calculate the ratio between two particles 
with the same chemical composition. The chemical fac- 
tors cancel out, and we are left with 



Y^ ^ fg2\ Sim^,pl) 
Y \ g J S{m^,p±_)' 



(18) 



where g* and g refer to each particle's degeneracy and 
the function S{m±,p±) is given by Eq. ((TBll . (Note that 
m± is the same for Y* and Y, but p± varies). 

Figure |3 shows the application of this procedure 
to the cases {K* + W)/{Ks) (top), S]*(1385)/A (mid- 
dle), and ri'/rj (bottom) at two freeze-out temperatures 
and flows: T = lAOMeV , Vmax / c = 0.55 on left and 
T — 170MeY , Vmax / c = 0.3 on the right. Significant de- 
viations from simple constant values are observed, show- 
ing the sensitivity of the ratio to freeze-out geometry and 
dynamics. The analytically simple result in Eq. H18|l is 
valid only if the light particle Y has been corrected for 
feed down from resonances, including Y* . In other words. 
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Eq. lO as well as Fig. |3 require that decay products 
from reconstructed Y* do not appear on the bottom of 
the ratio. Experiments usually do not do such feed down 
corrections !l9l since this would increase both 

statistical and systematic error on the ratio, and it is not 
always possible to do such corrections at all (undetected 
decays) or in the full range of experimental sensitivity. 

Introducing the feed down corrections into Eq. H18|l . 
we obtain 

Ypbserved ^ 9*S{mjL,pl) 

Yobserved gS{mj_,p±) + J^t diby -.Y Rim±, PTi) ' 

(19) 

Here, S{mj_,p±) describes the directly produced parti- 
cles and has the form given by Eq. (|15|l and each term 
R(rn±,p^^) describes a feed down contribution. 

In the case of an incoherent many-particle system, 
such as that we are dealing with, the dynamical (matrix 
element) part of the decay amplitude factors out [i^ . 
and i?(m^,p^,-) is obtained by integrating the statisti- 
cal hadronization distribution with a weight given by the 
phase space elements of the decay products. Thus, for a 
generic Y* ^ Y feed down is given by an iV-body decay, 

R{m^,p±) = (20) 

/ n ^S{m*r,PT)S (^p;-P, Ep,,^ , 

where the integral is performed over the whole allowed 
region. If more than one feed down occurs, Eq. (|2U|) can 
be used iteratively, with the left-hand side to be fed back 
to the right-hand side at each successive iteration. 

In general, this expression can get very complicated, 
and the Monte Carlo integration becomes necessary. For 
most cases considered here, where there is one feed down 
and two or three body decays, Eq. [501 can be carried out 
semianalytically [3 IM El • 

Figure 0| shows the ratios, including feed down of res- 
onances, for the same particles and statistical hadroniza- 
tion conditions as were studied in Fig.O In the S]*/(all 
A) case we omitted the feed down from S to A which is 
usually corrected for (if this is not done the ratio E*/ (all 
A) would depend strongly on the chemical potentials). 
We did allow for the (j) KsKl feed down, since it is a 
strong decay that cannot so easily be corrected for. We 
note that the feed down from particles with a different 
chemical composition cannot always be corrected for, and 
thus some resonances ratios will also acquire a (mild) de- 
pendence on the chemical potentials. This is even true 
for ratios such as r//(all ?]), given different ss content, in 
this paper, these type chemical corrections were set equal 
to unity. 

To further study the sensitivity of resonance-particle 
771^-ratio to freeze-out dynamics, we also present the 
(feed down corrected) case as a function of p± rather 
than TOj_ in Fig. [S] Unsurprisingly, we see grossly differ- 
ent behaviors, with many of the results coalescing. This 
of course is an expression of the fact that Y* and Y have 



dramatically different p±^ at the same m± and vice versa. 
We believe that the m± ratio will in general be more sen- 
sitive to freeze-out dynamics, since its dependence on m± 
is dominantly due to freeze-out geometry and dynamics. 
However, the p± dependence seen in Fig. |S1 provides an 
important self-consistency check for our previous results. 
We have found that the ■m± ratios are often greater than 
unity even though there must be more ground state par- 
ticles than resonances. Now it can be seen in the p± 
ratio, that this requirement is satisfied. 



IV. DISCUSSION 

In general the the mx and px dependence of the ratios 
in Fig 3 and, respectively, Fig 4 depends on freeze-out ge- 
ometry and dynamics. Changes in temperature and flow 
velocity alter the shape. The introduction of a steeper 
flow profile will further raise all of the considered ratios, 
since a considerable fraction of particles will be produced 
in regions that do not flow as much. The effect of freeze- 
out dynamics will generally go in the same direction as 
freeze-out approaches the explosive limit [dtf/dr 1). 
However, both the magnitude and the qualitative fea- 
tures of the two effects (flow and freeze-out velocity) will 
be considerably different. Especially, when more than 
one ratio is measured, it would appear that we will be 
able to determine the freeze-out condition. This is in 
contrast to the distributions in Fig. 13 where the ef- 
fects discussed in this paper result in linear corrections, 
which tend to compete, making the task of extracting the 
freeze-out dynamics much more ambiguous. Thus, there 
is considerable potential of resonance-particle mj^-ratios 
as a freeze-out probe. 

The presence of a long living hadronic gas rescatter- 
ing phase can distort our freeze-out probe. In particular, 
the apparent Y* /Y ratio will be altered due to the de- 
pletion of the detectable resonances through the rescat- 
tering of their decay products. It's dependence on mj_ 
will be affected in a non-trivial way, since faster (higher 
Pj^) resonances will have a greater chance to escape the 
fireball without decaying, thus avoiding the rescattering 
phase altogether. Regeneration of resonances in hadron 
scattering may add another m_L dependence which is dif- 
ferent for the E*/A and the K* /K ratios HO]. Other 
signals of the existence of such an interacting hadron gas 
phase have been considered 0,(23]. Fortunately, there is 
no evidence that a rescattering phase plays a great role 
in particle distributions. Even so, it would seem that 
the "safest" probes for freeze-out are the particles and 
resonances most unlikely to rescatter. 

For this reason we have included the rj' /rj ratio in 
our considerations, rj — > 77 and 77' — > 77 have very 
different branching ratios, but have the same degenera- 
cies and similar but rather small partial widths. The 
electromagnetic decay mode is practically insensitive to 
posthadronization dynamics. Regeneration effects are 
suppressed since the hadronic two body decay channel 
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is suppressed. All these features make these particles 
interesting probes, allowing for the analysis considered 
here. 77, ry' mesons have been measured at SPS energies 
in the 77 decay channel 0j \^ j and detectors such as 
PHENIX are capable of reconstructing the same decays 
at RHIC. 

While a long rescattering phase would affect the S* dis- 
tribution, the effect would be very easy to detect experi- 
mentally: 95% of S* decay through the p-wave S* — > Att 
channel. However, regenerating S* s in a gas of A s and 
TT s is considerably more difficoult, since Att scattering 
will be dominated by the s-wave Att ^ E^. This situa- 
tion will not occur for K* ^ Kit, since both decay and 
regeneration happen through the same process, leading 
to a very fast reequilibration time "so"]. Since both T,* /A 
and K* /K ratios have been calculated within the ther- 
mal model ^3 (neglecting rescattering), a strongly de- 
pleted S*/A ratio (compared with K* /K ) suggests that 
a statistical freeze-out description, such as that given in 
this paper, is incomplete, and an interacting hadron gas 
phase is also necessary. 



In summary, we have presented an overview of the dif- 
ferent statistical freeze-out models used to fit heavy ion 
data. We have shown how the freeze-out geometry and 
freeze-out dynamics influences the hadron spectra. Our 
primary result is the finding that the m± dependence of 
the resonance-particle ratios is a probe of freeze-out. We 
have presented these ratios for three particle species and 
two freeze-out conditions and have considered how our 
results could be altered by posthadronization phenom- 
ena. 



Acknowledgments 

Supported by a grant from the U.S. Department of En- 
ergy, DE-FG03-95ER40937. We thank Patricia Fachini 
(BNL) and Zhangbu Xu (Yale), from the STAR collab- 
oration, as well as Marcus Bleicher (ITP-Frankfurt) for 
fruitful discussions. 



[1] E. Fermi, Prog. Theor. Phys. 5, 570 (1950). 

[2] I. Pomeranchuk, Proc. USSR Academy of Sciences (in 
Russian) 43, 889 (1951). 

[3] L. D. Landau, Izv. Akad. Nauk SSSR (Ser. Fiz.) 17, 51 
(1953); English edition, L. D. Landau, Collected Papers 
of L. D. Landau, edited by D. Ter Haar, Pergamon, Ox- 
ford (1965). 

[4] R. Hagedorn, Suppl. Nuovo Cimento 3, 147 (1965). 

[5] J. Letessier and J. Rafelski, Hadrons and Quark-Gluon 

Plasma (Cambridge University Press, Cambridge, 2002). 
[6] P. Braun-Muntzinger, L Heppe and J. Stachel, Phys. 

Lett. B 465, 15 (1999). 
[7] L. Sandor, NA57 collaboration J. Phys. G in press 

, see http://wa97.web.cern.ch/WA97/Publications.html 

Hyperon production at the CERN SPS: results from the 

NA57 experiment^' 
[8] F. Becattini, M. Gazdzicki and J. SoUfrank, Nucl. Phys. 

A 638, 403 (1998); 

F. Becattini, J. Cleymans, A. Keranen, E. Suhonen and 

K. Redlich, Phys. Rev. C 64, 024901 (2001) 

F. Becattini, J. Phys. C 28, 1553 (2002). 
[9] S.V. Afanasiev et. al., NA49 Collaboration, Nucl. Phys. 

A715, 161 (2003). 
[10] G. Torrieri and J. Rafelski, New J. Phys. 3, 12 (2001) 
[11] S. V. Akkelin, P. Braun-Munzinger and Y. M. Sinyukov, 

Nucl. Phys. A 710, 439 (2002) 
[12] G. Torrieri and J. Rafelski, J. Phys. G 28, 1911 (2002) 
[13] J. Rafelski and J. Letessier, Nucl. Phys. A 715, 98, (2003) 

J. Letessier and J. Rafelski, Int. J. Mod. Phys. E 9, 107 

(2000) 

[14] J. M. Burward-Hoy, PHENIX collaboration, Nucl. Phys. 
A A715, 498, (2003). 

[15] J. Adams et al., STAR Collaboration, Multi-Strange 
Baryon Production in Au-Au collisions at ^/snn = 130 
GeV, nucbex/0307024 submitted to Phys. Rev. Lett. 

[16] K. A. Bugaev, M. Gazdzicki and M. I. Gorenstein, 



"hep-ph/0211337" 
[17] W. Broniowski and W. Florkowski, 

Phys. Rev. Lett. 87, 272302 (2001) 
[18] D. Magestro, J. Phys. G 28, 1745 (2002) 
[19] P. Fachini, STAR Collaboration, Nucl. Phys. A715, 

462 (2003); J. Phys. G in press, nucbex/0305034' 

Haibin Zhang, STAR collaboration, J. Phys. G in press, 

nucbex/0305034 
[20] G. Van Buren, STAR collaboration, Nucl. Phys. A A715, 

129, (2003). 

[21] C. Markert, STAR Collaboration, J. Phys. G 28, 1753 
(2002). 

[22] V. Friese, NA49 Collaboration, Nucl. Phys. A 698 (2002) 
487. 

[23] G. Torrieri and J. Rafelski, Phys. Lett. B 509, 239 (2001) 
and Phys. Rev. C64, 054907 (2001) Erratum-ibid. C65, 
069902 (2002). 
[24] E. Schnedermann, J. SoUfrank and U. W. Heinz, 

Phys. Rev. C 48, 2462 (1993) 
[25] E. Schnedermann, J. SoUfrank and U. W. Heinz, in Fire- 
ball spectra: NATO Advanced Study Inst.: Particle Pro- 
duction in Highly Excited Matter, II Ciocco, Italy, Jul 
12-24, 1992, H. Gutbrod and J. Rafelski, edts. proceed- 
ings series volume B303 (1993) 
[26] U. W. Heinz and P. F. Kolb, Proceedings of the 18th 
Winter Workshop on Nuclear Dynamics R. Bellwied, J. 
Harris, and W. Bauer, (EP Systema, Debrecen, Hungary, 
2002), pp. 205-216 hep-ph/0204061 
[27] F. Cooper and G. Frye, Phys. Rev. DIO 186 (1974) 
[28] K. A. Bugaev, Nucl. Phys. A 606, 559 (1996) 
[29] C. Anderlik et al, Phys. Rev. C 59, 3309 (1999) 
[30] D. Teaney, J. Lauret and E. V. Shuryak, Phys. Rev. Lett. 

86 4783 (2001). 
[31] T. Susa, NA49 Collaboration, Nucl. Phys. A 698 (2002) 
491. 

[32] M. Gazdzicki, NA49 Collaboration, proceedings of the 



8 



30th International Workshop on Gross Properties of [42] 
Nuclei and Nuclear Excitation: Hirschegg 2002: Ultra- 
relativistic Heavy Ion Collisions, Hirschegg, Austria, [43] 
13-19 Jan 2002 Edited by M. Buballa, W. Norenberg, [44] 
B.-J. Schaefer, J. Wambach. Darmstadt, Germany, GSI, 

2002. 377p. paper available as: 

Ihttp: / / theory.gsi.de /hirscliegg/2002 /Proceedings/ Gazdzickil 
137.138.139.124 draft l.ps.gz ' ' [45] 

[33] S. V. Afanasiev et ai, NA49 Collaboration, J. Phys. G 
28, 1761 (2002). 

[34] A. Olszewski et al, PHOBOS Collaboration J. Phys. G [46] 
28, 1801 (2002). 

[35] P. Staszel et al, Acta Phys. Polon. B 33, 1387 (2002). [47] 
[36] J. D. Bjorken, Phys. Rev. D 27, 140 (1983). 
[37] P. J. Siemens and J. O. Rasmussen, Phys. Rev. Lett. 42, [48] 
880 (1979). 

[38] J. Letessier and J. Rafelski, J. Phys. G 28, 183 (2002) [49] 
[39] J. Rafelski and J. Letessier, Phys. Rev. Lett. 85, 4695 
(2000) 

[40] E. Schnedermann and U. W. Heinz, Phys. Rev. C 47, 
1738 (1993). 

[41] S. Pratt, Nucl. Phys. A715, 389 (2003). [50] 
A. Dumitru and R. D. Pisarski, Nucl. Phys. A 698, 444 
(2002) 



A. Keranen, J. Manninen, L. P. Csernai and V. Magas, 
Phys. Rev. C 67, 034905 (2003) 

B. Touschek, Nuovo Cimento B 58, 295 (1968). 

R. Hagedorn, I. Montvay, and J. Rafelski, in: Hadronic 
Matter at Extreme Energy Density, N. Cabibbo and L. 
Sertorio, Edts. (Plenum Press, New York 1980), pp. 49- 
148. 

G. Torrieri and J. Rafelski, J. Phys. G in press and about 
to be published, nucl-th/0305071 Ullrich Heinz, private 
communication 

M. M. Aggarwal et al, WA98 Collaboration, Phys. Rev. 
Lett. 83, 926 (1999) 

R. Albrecht et al, WA80 Collaboration., Phys. Lett. B 
361, 14 (1995) 

A. Lebedev et al. [WA80 Collaboration.], Nucl. Phys. A 
566, 355C (1994). 

J. SoUfrank, P. Koch and U. W. Heinz, Prepared for 
HIRSCHEGG '91: 19th International Workshop On 
Gross Properties Of Nuclei And Nuclear Excitations. 
Edited by Hans Feldmeier. Darmstadt, Germany, GSI, 
1991. 243p. 

M. Bleicher, Nucl. Phys. A 715, 85, (2003). 



9 




T=140 MeV 
v=0.55 



no flow 

dt/dr=1 , no spread 
dt/dr=1 , y~r 
dt/dr=1 , \i~r' 
Hubble, no spread 
Hubble, spread {iT-m^) 
dt/dr=0, no s 
dVdr=0, y~i 
dt/dr=0, v~i 



4.0 



3.0 



2.0 - 



1.0 - 



T=170 MeV 
v=0.3 



no flow 

dt/dr=1 , no spread 
dt/dr=1 , y~r 
• dt/dr=1 , v~r' 
Hubble, no spread 

■ Hubble, spread (L^'- 
dt/dr=0, no spread 
dt/dr=0, y~r 

■ dt/dr=0, v~r' 



1.0 
6.0 



1.5 



2.0 
[GeV] 



2.5 



3.0 



0.0 



1.0 



1.5 



2.0 
[GeV] 



2.5 



3.0 



4.0 



2.0 



0.0 




T=140 MeV 
v=0.55 



1.3 



1.8 



2.3 
[GeV] 



no flow 

dt/dr=1 , no spread 
dt,/dr=1 , y~r 
dt,/dr=1 , v-r" 
Hubble, no spread 
Hubble, spread 
dt/dr=0, no spread 
dt,/dr=0, y~r 
dt,/dr=0, v~r' 



2.8 



6.0 



4.0 



T=170 MeV 
v=0.3 



no flow 

dt/dr=1, no spread 
dt/dr=1, y~r 
■ dt/dr=1, «~r' 



Hubble, no spread 

Hubble, spread (e'-tj'') 

dt/dr^O, no spread 

dt/dr=0, y~r 

dt/dr=0, «~r' 



2.0 



1.3 



1.8 



2.3 
[GeV] 



2.8 





T=140 
v=0.55 



no flow 

dt/dr=1 , no spread 
dVdr=1, y~r 

dt/dr=1,v~r' 
Hubble, no spread 
Hubble, spread 
dt/dr=0, no spread 
dt/dr=0, y~r 
dt/dr=0, v-r" 



0.5 



0.0 



1.0 



1.5 



2.0 
[GeV] 



2.5 



3.0 



0.0 



T=170 MeV 
v=0.3 



no flow 

df/dr^1, no spread 
dt/dr=1, y-r 

• dt/dr=1, «~r" 
Hubble, no spread 

■ Hubble, spread if-v,^] 
dt/dr=0, no spread 
dt/dr=0, y~r 

• dt/dr=0, «~r' 



1.0 



1.5 



2.0 
[GeV] 



2.5 



3.0 



FIG. 3: (Color online) Dependence of the K* /K, E*/A and r;' /eta on the Freeze-out model 
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FIG. 4: (Color online) (K* + K )/(all Ks), E* (1385)/(all A) and ?7'/(all r\) ratios, including feed down from resonances. 
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FIG. 5: (Color online) dependence of (K* + K )/(all Ks), E* (1385)/(all A) and rj' /{all rj) ratios, including feed down from 
resonances. 



